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ABSTRACT 

In order to study the thermodynamic properties of brane-antibrane systems 
in the toroidal background, we compute the finite temperature effective potential 
of tachyon T in this system on the basis of boundary string field theory. We 
first consider the case that all the radii of the target space torus are about the 
string scale. If the Dp-Dp pair is extended in all the non-compact directions, 
the sign of the coefficient of |T| 2 term of the potential changes slightly below 
the Hagedorn temperature. This means that a phase transition occurs near the 
Hagedorn temperature. On the other hand, if the Dp-Dp pair is not extended in 
all the non-compact directions, the coefficient is kept negative, and thus a phase 
transition does not occur. Secondly, we consider the case that some of the radii 
of the target space torus are much larger than the string scale and investigate 
the behavior of the potential for each value of the radii and the total energy. If 
the Dp-Dp pair is extended in all the non-compact directions, a phase transition 
occurs for large enough total energy. 
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1 Introduction 



Non BPS configurations of branes such as brane-antibrane pairs and non BPS D-branes have 
received attention recently (for a review see e.g. [1]). Open strings on these branes have 
a tachyonic mode, and these configurations are unstable. The tachyon potential on these 
configurations is calculated by using string field theory. For example, the tachyon potential 
on the Dp-Dp system in type II string theory was computed on the basis of boundary string 
field theory (BSFT) [2] [3] by Minahan and Zwiebach [4] and by Kutasov, Marino and Moore 
[5]. It is given by 



where T is a complex scalar tachyon field, V is the volume of the system that we are consid- 
ering, and r p is the tension of a Dp-brane, which is defined by 



where g s is the coupling constant of strings. This potential has the minimum at |T| = oo, 
and its qualitative features agree with Sen's conjecture [6]. 

Although the tachyon potential is unstable at T = 0, there is a possibility that 
T = becomes a stable minimum at high temperature. There are some approaches to 
the statistical mechanics of open strings on the Dp-Dp system. Danielsson, Giiijosa and 
Kruczenski evaluated the finite temperature effective potential when there is only the tachyon 
field [7]. Huang attempted to evaluate the finite temperature effective potential based on 
BSFT of bosonic strings [8] and then investigated a finite temperature system of the parallel 
brane and antibrane with a finite distance [9] . However, in these models, we cannot see the 
behavior of the system near the Hagedorn temperature, which is the maximum temperature 
of perturbative strings. 

In the previous work [10], we have computed the finite temperature effective potential 
of the brane-antibrane systems near the Hagedorn temperature on the basis of the work of 
Andreev and Oft [11]. In the D9-D9 case, the sign of the coefficient of |T| 2 term of the 
potential changes slightly below the Hagedorn temperature from negative to positive as the 
temperature increases, while the coefficient is kept negative in the Dp-Dp case with p < 8. 
Therefore, we conclude that not a lower dimensional brane-antibrane pair but a spacetime- 
filling D9-D9 pair is created near the Hagedorn temperature. This is the result in the case 
that the target space is non-compact. 

A considerable number of studies have been made on the application of the brane- 
antibrane system to cosmology [12] [13] (for a review see e.g. [14]). It is important to 
investigate the brane-antibrane systems from the cosmological viewpoint because the tem- 
perature is extremely high in the early universe, and brane-antibrane pairs might exist in 
this period. Then the universe expands inflationary, because the tension energy of the brane- 
antibrane pairs can provide an effective cosmological constant. In the early universe, the size 
of the universe is extremely small and the topology of the universe is also very important. 
This is because Strings and branes can wrap the compact directions of target space and they 
are sensitive to its topology. The purpose of this paper is to investigate the brane-antibrane 
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system at finite temperature in the toroidal background, which is the simplest non-trivial 
one. 

This paper is organized as follows. In §2 we evaluate the free energy of open strings 
on the Dp-Dp system in the toroidal background on the basis of BSFT. We must compute the 
finite temperature effective potential by using the microcanonical ensemble method, because 
we cannot trust the canonical ensemble method near the Hagedorn temperature as we will 
explain later. Then we investigate the behavior of the free energy in the complex /3-plane 
in §3 and compute the finite temperature effective potential in the microcanonical ensemble 
method in §4. §5 presents our conclusions and discussions. 



2 Free Energy of Open Strings 

It is well-known that we must compute the statistical variables of strings near the Hagedorn 
temperature by using the microcanonical ensemble method [15]. This is because we cannot 
trust the canonical ensemble method near the Hagedorn temperature from the following 
reason. The partition function Z(/3) is given by the Laplace transformation of the density 
of states Q(E): 

poo 

Z(f3) = / dE Q(E)e~ /3E . (2.1) 
Jo 

We can expect that the canonical ensemble method gives the same quantities with that of 
the microcanonical ensemble method if the integrand has a sharp peak. However, the density 
of states of strings behave as 

Sl(E) ~ e feB , (2.2) 

for large E, where f3 H is the inverse of the Hagedorn temperature, and the integrand has no 
sharp peak near the Hagedorn temperature. Thus, we must compute the finite temperature 
effective potential by using the microcanonical ensemble method. Q(E) can be obtained 
from the inverse Laplace transformation of Z(/3): 

Q(E) = [ L+10 ° ^Z([3y E . (2.3) 

The density of states has been derived from this formula in the case of closed strings [16] 
[17] and in the case of open strings on D-branes [18], and in the case of open strings on the 
Dp-Dp system [10]. We will apply this method in the case of open strings on the Dp-Dp 
system in the toroidal background. Since Z(f3) is related to the free energy F(f3) as 

Z(/3) = exp[-/3F(/3)], (2.4) 

we will begin by evaluating the free energy of open strings. 

Let us make a comment here. One may wonder why we use such an indirect method 
to compute the density of states. The density of states is defined by 

Q(E) — tr8(E — H) , (2.5) 
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where E is the total energy and H. is the Hamiltonian. It is difficult to take the trace over 
the states of strings directly. If we use complex integral representation of the delta function 
and the definition of Z(f3), which is given by 

Z((3) = tre~ m , (2.6) 

we obtain the inverse Laplace transformation (2.3) [19]. 



In order to compute the free energy by using the Matsubara method in the ideal gas 
approximation, we must compute the one-loop amplitude of strings. Much work has been 
done to generalize BSFT to the one-loop level of open strings on the Dp-Dp system [20]. 
But there is an ambiguity in choosing the Weyl factors of the two boundaries of a one-loop 
world-sheet, because the conformal invariance is broken by the boundary terms in BSFT. 
Andreev and Oft have proposed the following form of 1-loop amplitude [11] in type II string 
theory from the principle that its low energy part should coincide with that of the tachyon 
field model [4] [5]. If we restrict ourselves to the constant tachyon field, which we denote by 
T, it is given by 



IQttHV f°°dT /A ._„+i 
Z t = - -^+r / — 4vrr 2 



-A-k\T\ 2 t 



(2tto 



X 



A'(o|«-), 



' fl 2 (0|zr) ' 
v 0i'(O|ir). 



(2.7) 



where V is the volume of the system that we are considering. Here we are considering only 
the coincident brane-antibrane system, and we will not treat the parallel brane and antibrane 
with a finite distance [21] [9] in this paper. We can obtain the same amplitude when we 
consider the field theory which has the mass spectra 



M NS 2 = -(N B + N NS + 2\T\ 2 
a' v 



2) ' 



Mr 2 



}- i (n b + N r + 2\T\ 2 ), 



(2.8) 
(2.9) 



where M NS and M R are the mass of the Neveu-Schwarz and Ramond sectors, respectively, 
and N B , Nns and Nr are the oscillation modes of the boson, Neveu-Schwarz fermion and 
Ramond fermion, respectively. If we use the proper time form of free energy [22], which is 
given by 



= a£ ±i L — ( 47rT ) 2 

{2-na')— Jo t 



]T X>xp \-2ita'M NS 2 T - 



M NS 



2 r= l 



4sr r d i^y v ^ E E(-i) r - P f-2W 

■a') 2 Jo r =1 V 



Mr 2 t- 




87ia'T 



(2.10) 



we can compute the free energy. The result is 

16tt 4 V r 00 dr 



F(T,(3) 
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0i'(O|ir) 
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tf 2 (0|fr) 



This is the free energy on the non-compact background [10]. 
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Let us generalize this free energy to that in the toroidal background. We consider 
the Dp-Dp system in type II string theory compactified on D-dimensional torus Tp. In the 
compact direction, there are two kinds of zero modes which give discrete energy spectra. 
They are momentum modes in the directions parallel to the Dp-Dp system, and winding 
modes in the directions transverse to it. If we include these modes to the mass spectra, it is 
given by 



M NS 2 



Mr 2 



p—d 
p—d 



D 



=p-d+l 



a' 



a' 



+ -[N B + N NS + 2\T\ 2 - - , (2.12) 



l=i v-n.// i= p -d+i v « / « 



Here, we assumed that the Dp-Dp system is extended in the <i-dimensional non-compact 
directions. It should be noted that these mass spectra are invariant under T-duality trans- 
formation [23] 



Ri 



a' 
Ri 



for the directions parallel to the Dp-Dp system and 

j 



Ri 



a 
Ri 



Hi 



(2.14) 



(2.15) 



for the directions transverse to it. Summations over the momentum modes and winding 
modes can be rewritten by using i? function as 



E 



exp 



mj,rii=— oo 



p—d , * 9 



D 

+ E 

i=p— d+1 



riiRj 



a' 



p-d D / 

= n n 4 

7=1 i= p -d+l \ 



^4 



2iR 2 r ' 



Thus, the free energy is given by 
F(T,P,R) = 



f^f i-^ff n 4 

Ph Jo t 2 i=n= P -d+i \ 



a' 



2ia'r 



(2.16) 



Ri 2 



A 



2%R 2 t 



a' 



x 



^i'(0|ir) 



t? 3 



iP 2 



' MOW 

MOW. 



- 1 



i/3 2 



(2.17) 
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where Vd is the <i-dimensional volume in the non-compact directions parallel to the Dp-Dp 
system, and f3 H is the inverse of the Hagedorn temperature 

P„ = 27r v / 2a 7 . (2.18) 

This free energy is invariant under the T-duality transformation. The finite temperature 
effective potential is also invariant under the T-duality transformation, because it is derived 
from this free energy. We only need to investigate the region R > y/a', and we restrict 
ourselves to this region hereafter. In the next section, we will investigate the singular part 
of the free energy in the complex /3-plane. 



3 Complex Temperature Formalism 



Let us investigate the behavior of the free energy near singularities in the complex /3-plane. A 
general property of the free energy of strings is that it has the leading singularity at f3 — f3 H , 
which we call the Hagedorn singularity. If some radii are much larger than the string scale 
\[a! , there is the case that other radius dependent singularities exist close to the Hagedorn 
singularity [16] [17] [18]. We will show these properties of the singular part of the free energy 
in this section. 

In order to see the behavior of the free energy at high temperature, let us make the 
variable transformation 

1 

T = ? 

in (2.17). Using the modular transformation of -Q functions, we obtain 



F(T,P,R) = 



p-d p 2 



2-& 

f 
Jo 



|2\ P-d D 



f°° , D+d-9 I All IT, , „ , 

x / dt ts- exp ( L-L ) J] J] ^ 3 



x 



v 1 

mm) rl° 
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1? 
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2a' 



#3 



ia't 
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^4 



- 1 
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(3.1) 



Expanding d functions and extracting the leading term in the large t region near the Hage- 
dorn singularity, we obtain 



F(T,P,R) 



D+d-9 



dt t 2 exp 



p-d D / 

xn n mo 

i=i i =p -d+i \ 
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A 
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(3.2) 
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where we have defined 



and 



A = 



a 



,d-p+% 



p—d 

v P = v d n Ri, 

1=1 



(3.3) 



(3.4) 



and we have introduced the low energy cutoff A. This cutoff is required because we cannot 
apply above approximation to the low r region. Since we will compute the |T| 2 term of the 
finite temperature effective potential in the vicinity of T = 0, let us expand the free energy 
in |T| 2 and keep the lower order terms 



F(T,(3,R) ~ 



dt t 



D+d-a 



Ph Ja 
p-d D I 

*n n m 

i=i i =P -d+i \ 

47v\T\ 2 AV n r° 



exp — 7T- 







%Ri 2 t 



P 2 ~P^_ 

2 
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Pn 



t 



2a' 
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2R> 



1 H Ja 



P. 

p-d D / 
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D+d-ll 

at t 2 exp 



-7T- 



P 2 ~ Ph 



p, 



iRjH 



2a' 



ia't 



2Ri 2 



(3.5) 



We can evaluate the singular part of the free energy by computing the discontinuity around 
the singularity like in the previous paper [10]. But let us compute it by using the incomplete 
r function for the later convenience. 



3.1 String Scale Radius Case 

If all the radii are close to the string scale \fa\ the Hagedorn singularity is dominant and 
we only need to consider the terms 

+ 4n\T\ 2 AV p [°° J± ^ D+d-u ___ ( _p 2 - p 2 



'-Tdtt^exp ~n^4^t . (3.6) 
Ja \ Ph J 



P H Ja "V P 

This free energy is proportional to that in the non-compact background case [10] if we replace 
D + d by p. From this we can expect that a phase transition occurs only when D + d = 9 
and in fact it is true as we will see in the next section. We can rewrite (3.6) near P — P H by 
using the incomplete T function as 



n 



2(2^AV P (P - p H )^\T ? r / ^ 2 J_-p^ k \ ^ 



P H a+2 V Ph 
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where we have defined 

7 -D-d 

a = . (3.8) 

We can deduce the singular part of the free energy F sing from (3.7) as follows. 



(a) D + d = 9 (a = -1) 

When D + d = 9, which means a = —1, the first argument of the incomplete T function 
in the first term of (3.7) becomes one, so that we can set A = 0. For the second term 
of (3.7) we can use the following formula for the incomplete T function; 

00 (—l) n x n 

T(0, x) = - 7 - log* - £ 1 > (3.9) 

n=l n ' n - 

where 7 is the Euler constant. Combining these two terms, we get 

We can take the same notation as the previous work [10] by setting A = (2tt)~ 1 . We 
will set A = (27r) _1 in this formula hereafter. 



(b) D + d : even (a : half- integer) 

When D + d is even, the first arguments of the incomplete T functions are negative 
half-integers. In this case we can set A = and we obtain the singular part of the free 
energy as 



2(*)-«IX-a-lMl 4 jTP 



a+2 



(c) D + d : odd (a : integer) 

When D + d is odd, the first arguments of the incomplete T functions are negative 
integers. The incomplete T function whose first argument is a negative integer can be 
expanded as 

1 n (—l) n 
T(-n, x) = —e~ x ^(-l)*- 1 ^ - s)\ x-^ 8 - 1 + ^— p-r(0, x). (3.12) 

In our case, A — > corresponds to x — > 0, so that the incomplete T function can 
be approximated by the last term. Thus, the singular part of the free energy can be 
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Figure 1: Radius dependent singularities 



obtained as 

F smg (T,(3,R) * - r{a + 1)(3H a + i (P-P*) log ^2*-^— A 



.l)-2(2.)-AV P |T| 2 (/5 _^ r+llo / 2 ^ A ^ 



r(« + 2)/3 H Q+2 "" ~°v~" & 

(3.13) 

We will also set A = (27r) _1 in this formula hereafter. 



3.2 Large Radius Case 

Let us consider the case that some radii are much larger than the string scale \fa'. When the 
radii in the directions parallel to the Dp-Dp system are much larger than the string scale, 
the Hagedorn singularity is dominant and we can use the results above. But when the radii 
in the directions transverse to the Dp-Dp system are much larger than the string scale, the 
radius dependent singularities exist near the Hagedorn singularity as is sketched in Figure 
1, and we must consider the radius dependent terms. We will show only the D + d = 9 case, 
because, in the D + d < 8 case, the result can be obtained in a similar way to the D + d = 9 
case as we will comment in the next section. For simplicity, let us assume that some of the 
transverse directions has the same radius Rq with Rq 3> \J~ol' . Since ■d function in (3.5) can 
be rewritten as 



p—d+d / 

n M° 

i=p-d+l \ 



ia't 



9R _2 j -= L ex P " L -^Dr*l> ( 3 - 14 ) 

zn -d I n;=-oo \ i=p-d+l ZrL d 



the free energy can be approximated as 



F(T,P,R) ~ --p/ dt 2 ex P 



P 2 - p n 2 

-7T 




H s /. T 2^ ex P _7r — t ' ( 3 - 15 ) 



where we have defined 



p—d+d 2 I 

ma 



Pn 2 = P H 2 \l- E (3-16) 

i=p-d+l An d 
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From this we can see that singularities exist at (3 — /3 n on the real axis of complex /3-plane, 
and they approach to the Hagedorn singularity (3 = (3 H as Rj increases. This free energy can 
be rewritten by using the incomplete Y function as 



F(T,(3,R) ~ £ 



n^=— oo L 



(3 H AV P 



AnAV p \T 



- r (i , 2^; a ^A 

- r f , 2«A 



(3.17) 



where we assumed that (3 ~ /3 n for each terms. For the first term, we can set A = 0. For the 
second term, we can use (3.9). Combining these two terms, we get 



F sing (T,(3,R)~ £ 



B H AV P 



4irAV p \T 



2n(3 n ((3 - (3 n ) 



(3 L 



P, 



(3.18) 



We will also set A = (2n) in this formula hereafter 



Finite Temperature Effective Potential near the Hage- 
dorn Temperature 



In this section, we compute the finite temperature effective potential of the brane-antibrane 
systems in the toroidal background. We must compute it by using the microcanonical en- 
semble method as we have mentioned in §2. All the statistical variables are derived from 
the density of states in the microcanonical ensemble method. The density of states Q(E) is 
obtained from the inverse Laplace transformation of the partition function Z((3) (2.4). The 
partition function is given by (2.4), which is rewritten as 

Z{(3) = exp [-(3{F reg + F sing )\ , (4.1) 

where F reg is the regular part of the free energy. The entropy S(E) can be derived from 
Q(E) as 

S(E) = log Q(E)5E, (4.2) 

where 5E represents the energy fluctuation, and the inverse temperature (3 from the partial 
derivative of S(E) 

*-£■ _ 

The finite temperature effective potential of the Dp-Dp system can be calculated from these 
variables as 

V(T,E) = 2r p V p exp(-8\T\ 2 )-p- l S. (4.4) 

Now, we are ready to calculate the finite temperature effective potential by using the micro- 
canonical ensemble method. 
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4.1 String Scale Radius Case 



Let us first consider the case that all the radii in the transverse directions are about string 
scale. In this case, the calculation of the finite temperature effective potential in the toroidal 
background is similar to that on non-compact background. Thus, we will calculate in the 
D + d = 9 and the D + d — 8 cases explicitly, and show only the results for the other cases. 



(a) D + d = 9 (a = -l) 



Let us consider the D + d = 9 case, which is the most interesting case. This is the case 
that the Dp-Dp system is extended in all the non-compact directions. In this case, 
F S ing is given by (3.10) and F reg can be expanded as 

-(5F reg {(3) = A V P - a V p (P - (3 H ) + 0{V p (P - (3 H ) 2 ). (4.5) 

Then, by using (3.10) the inverse Laplace transformation (2.3) can be rewritten as 



n(T, E, R) ~ e ^ E+x » v * [ ( ) ' exp 

Jc a 2m \ p H ) 



£>_iV p |T| 



2 



(P-P H )E + 



p 



p-p t 



(4.6) 



where we have defined 



D_, = -AttA, (4.8) 

and we take the contour C a as sketched in Figure 2. Let us suppose that both E and 
E/Vp are very large. Then, the saddle point method works well because the exponent 
in the integrand is very large. The result is 



I / 3 2 E \ 5 d - iV pI t I" ( r< y\l / , — s 

n(T,E,R)~ — l±^-y) exp [ PhE + XoVp + 2 v c -^p e ) ■ 

(4.9) 



Substituting Q(E) into (4.2), we obtain 



f-yv) 4 \C_JVH5E)I, 

+P H E + \ V P + 2^/c^VpE, (4.10) 

and from (4.3), we get 



H 2E AE 1 V E K J 

The finite temperature effective potential can be derived from (4.4). In order to argue 
the stability of the brane-antibrane system, we need only \T\ 2 term of V{T, E, R). This 
term is given by 



27rAV Pl (2ixp H E 



-16t p V p + —-^ log 



(4.12) 
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Figure 2: Complex /3-plane 



It should be noted that the second term in the coefficient of |T| 2 increases with in- 
creasing E. Since the first term is constant as far as V p and r p fixed, the sign of the 
|T| 2 term changes from negative to positive at large E. The coefficient vanishes when 

- AV P (8f3 H r p \ 

E ^2^- H exv [-^r)- (4 - 13) 

Since (4.11) can be approximated as 



(3^{3 H + J^, (4.14) 



H 



1 + exp 



(4.15) 



for large E, we can derive the critical temperature T c at which the coefficient vanishes. 
The result is 

irA 

Here we see that this temperature is very close to the Hagedorn temperature since 
t p is very large if the coupling of strings is very small. Above this temperature, the 
coefficient of \T\ 2 is positive and T = becomes the potential minimum. This implies 
that a phase transition occurs at the temperature T c which is slightly below the Hage- 
dorn temperature, and the Dp-Dp system, which is extended in all the non-compact 
directions, is stable above this temperature. 

(b) D + d = 8 (a = -1/2) 

Let us consider the D + d = 8 case, where we will see that a result strikingly different 
from the D + d = 9 case arises. This is the case that the Dp-Dp system is extended 



12 



in all the non-compact directions except for one direction. In this case, the inverse 
Laplace transformation (2.3) is rewritten as 



fl{T,E,R) ~ e PHE+x v P 



dp 
c a 2m 



cxp 



-DiV p \T\ 2 ([3-[3 H f 



~ e 



J3„E+\ V P 



dp r 
c a 2ni 

x exp 



l-D _MT\ 2 {P-p H )* 



(P-P H )E + CiV p (P-p H )-* ,(4.16) 



where we take the small |T| 2 approximation in the second equality and we have defined 



C i 

~2 



D i 



Pi 



(4.17) 
(4.18) 



We can also use the saddle point method and obtain 



Sl(T,E,R)~ ~\ ^ exp 

32237T2E« 



3C_izVp 3 E* 
P H E + A V P + - 1 



V 



2§ 



C_i3D_iV P 3 |T| 
g g f_ 1 1 



(4.19) 

We can calculate the entropy 5, the inverse temperature (3 and the potential V(T, E, R) 
from (4.2), (4.3) and (4.4) as in the D + d = 9 case. The |T| 2 term of V(T,E,R) is 
given by 



-16r p V p - 



2 5 7T 2 A3V/ 



\T\ 



(4.20) 



It should be noted that the second term in the coefficient of |T| 2 decreases as E gets 
large. Thus, the coefficient of |T| 2 remains negative for large E. This implies that a 
phase transition does not occur unlike in the D + d = 9 case. 



(c) D + d < 7 {a < 0) 

We only show the results for the D + d < 7 case, because we can obtain these results by 
the similar calculation to the previous work [10]. The |T| 2 term of the finite temperature 
effective potential is given by 



8tt 2 A 2 V 2 f/^F 
ToTpVp _ ^ log 



p / j 



I I ^ 

I I 5 



(4.21) 



for D + d = 7, 



2 10 ttM 4 V 4 
-16r p V p + —3^ 



(4.22) 
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for D + d = 6, 

T| 2 , (4.23) 
for D + d = 5, where 

O77-4 

Ef = E-— V p log(AV p ), (4.24) 

and 

T| 2 , (4.25) 

for .D + d < 4. From these we can see that the coefficients are negative for large E, so 
that a phase transition does not occur. 



4.2 Large Radius Case 

Next, let us consider the case that some of the radii in the directions transverse to the 
Dp-Dp system are much larger than the string scale. In this case, the radius dependent 
singularities exist near the Hagedorn singularity. The behavior of the finite temperature 
effective potential is determined by the relation between the total energy and the radii of 
the target space torus. We will mainly consider the D + d = 9 case. We assume that all the 
(i-dimensional directions have the same radius Rj, for simplicity. 



16tt 2 

- 16TpVp + p7e> 



-16t p V p + 



16tt 2 

~p7b 



(a) Hagedorn Singularity Dominance, Saddle Point Method 



In order to estimate the contribution from the radius dependent singularities, let us 
pick up the leading two singular terms of the free energy (3.18). Then the inverse 
Laplace transformation (2.3) can be rewritten as 



tt(T,E,R) ~ / — exp 



dp_ 

c a 2m 



(4.26) 



where C_ x is defined by (4.7) and p\ is the nearest singularity to the Hagedorn singu- 
larity, namely, f3 n with J2i n i 2 = 1- In this case, we can deform the contour as sketched 
in Figure 3, and we obtain 



Cf3 H 

d(3 



(3E - (3F reg {(3) + 



c.y v 



+ / 7T— exp 



PE-[3F reg (f3) + 2d 



P-Pl 

Let us expand F reg (f3) in a power of (f3 — f3 H ) as (4.5) and of {(3 — Pi) as 
-PF reg (P) = A X V P - a{V p {P - Pi) + 0(V p (P - px) 2 ), 



(4.27) 



(4.28) 
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Figure 3: Contour deformation 



and define 



E 



E - cxiVp, 



E 



(P-P 



HJ, 



z = 



E 



\ 2dc_y p 



(P-Pi 



Then (4.27) can be rewritten as 
n(T,E,R) ~ 



/3„E+X Vp 



E 



dz 
2tt7 



exp 



c_y p E [z + ~ 



2dC_ 1 V p c f3 lE+XlVv f d^_ 
\ E' J 2m 



exp 



2dC_ 1 V v E z' + 



(4.29) 



(4.30) 



where the contours for each terms encircle the origin in the complex z-plane and z'- 



plane counter-clockwise respectively. If we assume that yC^VpE ^> 1, then we can 
use the saddle point method. Substituting (4.7), we can rewrite this condition as 



By using the saddle point method, we obtain 

Q(T, E, R) ~ Q (T, E, R) + fii(T, E, R), 



(4.31) 



(4.32) 
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where Qq(T, E, R) and Qi(T, E, R) represent the contribution from the Hagedorn sin- 
gularity P = (3 H and from the next leading singularity (5 — Pi respectively: 



; l r< v \ 4 



n (T,E,R) ~ ^(^pLj exp{p H E + \ V p + 2yJC_ 1 V p E), 
_i i 

Qi(T, E, R) ~ f^r^ 4 exp + AxV p + 2^260^$} . 

2 4 V^ \ E J ^ J 

If Iog(f2 /f2i) ^> 1, then Q(T,E,R) is dominated by Q (T, E, R) and we can ignore 
the term Q 1 (T,E,R) [16] [18]. Let us assume that E > a V p . Then log(Q /^i) can 
be approximated as 

log^ ~ [(5 H E + A V P + 2^C^V P E^ - (prf + XiV p + 2^2dC_ 1 V p E') 

~ (P H - Pi)E - (Jid - l) 2ylc_y p E, (4.33) 

where, in the second equality, we have used (Ao — Ai)V p — — <ToV p (P H — Pi), which is 
derived from (4.5) and (4.28). Then above condition is approximated as 



{P H -P l )E^>^C_ 1 V p E. (4.34) 
Substituting (4.7) and (3.16) with Y^i^i 2 = 1, we obtain 

p_4-<Z-\ ) 

E^> 1 - p . (4.35) 
a' 2 

We can also show that we can ignore the contribution from other singularities. There- 
fore, if E satisfies 



E > max I a V p , * d * d Vl 



V? ' 5+p-d 



(4.36) 



then the density of states is dominated by the contribution from the Hagedorn singu- 
larity and we can use the saddle point method to compute the density of states. By 
using the saddle point method, we can also show the Hagedorn singularity dominance 
in |T| 2 terms. As a consequence, V(T, E, R) is the same as that in the case of d = 
and D + d = 9, and phase transition occurs near the Hagedorn temperature. 

(b) Hagedorn Singularity Dominance, No Saddle Point Method 

Let us consider the region where we cannot use the saddle point method. From (4.30) 



the saddle point approximation does not work well when JC^VpE <C 1, namely, 



p—d—l -j 

E « Ud . (4.37) 
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In this case, the exponents in the integrands in (4.30) are very small and both integrals 



are given by powers of yC_ 1 V p E. Thus, log(Q /^i) can be approximated as 

log^-O^-ft)^, (4.38) 

we can ignore the term £li(T,E,R) if E satisfies (f3 H — Pi)E 3> 1, which is rewritten 

as 

D_2 

E > -4-. (4.39) 
a' 2 

We can also show that we can ignore the contribution from other singularities in this 
region. Therefore, if E satisfies 

^oV p , « E « a ' ' Rl , (4.40) 

then the density of states is dominated by the contribution from the Hagedorn singu- 
larity and we can ignore the contribution from other singularities. Thus, the density of 
states is given by (4.6), which can be rewritten by using the modified Bessel function 
of the first kind as [18] 

n(T, E, R) ^ 2 p H »-iV P \Tf e p H E + x v P lD _ iVp ^_ 1 ( 2V / C r ^l) . 

(4.41) 

Let us expand the modified Bessel function as 

W-g o n!r(i/ + n+1) y • (4-42) 



and keep the lower order terms because the argument 2yC_ 1 V p E <C 1 is very small. 
Using the fact that the T function can be approximated as 



r(e) ~ 1 (4.43) 



for small e, we obtain 



|2 X 



n(T, E, R) ~ C_ 1 V p {(3 H E) D - lV - m2 e^ E+XoV - (l + ^= j • (4-44) 

We can calculate the entropy S, the inverse temperature f3 and the potential V(T, E, R) 
from (4.2), (4.3) and (4.4) as in the string scale radius case. The |T| 2 term of V(T, E, R) 
is given by 

\T\ 2 . (4.45) 



Air AV 

-16r p V p + — -^\ g(p H E) 

Ph 



The coefficient vanishes when 



^-L eX p[ 2 2 Ph + - TpR «' 1 , (4.46) 



Tea' 2 
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where we use the explicit expression of A with D + d — 9, that is, 



d-\-d — p 



A = 



a' 



. 5-d 



(4.47) 



But this E cannot satisfy the condition (4.37). This is because r p is very large if 

p , — 
the coupling of strings is very small, and V p > a 2 since we assumed that R > y a' '. 

This means that, as E increases, E exceeds the region (4.40) before a phase transition 

occurs. Thus, a phase transition does not occur in this region. 

(c) Contribution from Radius Dependent Singularities 

In the previous two regions (4.36) and (4.40), the density of states is dominated by the 
contribution from the Hagedorn singularity. Here, we consider the remaining region 



or 




p—d—1 "T ^ 



max 



o"oV p , 



^'<£< d p 



V 



5+p-d 



(4.48) 



(4.49) 



a' 2 



in which we cannot ignore the contribution from the radius dependent singularities. 
The singular part of the free energy F sing is given by (3.18). Although we must perform 
the summation over rii in order to include the contribution from the radius dependent 
singularities, it is difficult to compute this sum directly. However if we return to (3.15), 
which is rewritten as 



F(T,P,R) 



Ay roo 00 

8 L dt £ 

Pa JA „ l= _ c 



cxp 



p—d+d 2 I 

7m,- a 



E 



i= p -d+i ^-Rd 



t — 7T n 1 



+ 



4ir\T\ 2 AV v r°°dt 



Pe 



" a 7 ^ exp 



- E 



d+d imM 



, p-d+ I. 



(3 2 -(3 H . 

t — 7T — n 1 



(4.50) 



the sum over rii is approximated by integral as 



p—d+d oo 

n e ex p 

i=p— d+1 «i = -oo 



2ifc 



p—d+d 

n 



ixa't / j-oo 



(4.51) 



d / i=p—d+\ 

for large Rj, and we can compute it. By evaluating the gaussian integral, we obtain 



F(T,(3,R) ~ J~ dt H (^?-^t 



An\T\ 2 A'V p - 



d+1 



dt t 2 exp —n 



(4.52) 
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where we have defined 



a' p+ 2 



A = = . 4.53 

2^^S d u 9-d d. 

This free energy is equivalent to (3.6) with D + d = 9 — d. If we compute the finite 
temperature effective potential from this, we obtain the same results with the string 
scale radius case with D + d = 9 — d. As a consequence, a phase transition does not 
occur in this region. 



Here, let us make a short comment on the D + d < 8 case. For extremely large E, we 
can show that the contribution from the Hagedorn singularity is dominant, and we obtain 
the same results with the string scale radius case with the same D + d. For not so large 
E, the summation over rii can be approximated by the gaussian integral, and we obtain the 
same results with the string scale radius case with D + d replaced by D + d — d. Therefore, 
a phase transition does not occur in this case. 



5 Conclusion and Discussion 



In this paper, we have discussed the behavior of the finite temperature effective potential 
on the brane-antibrane pair in the toroidal background. We have evaluated the |T| 2 term of 
the potential near the Hagedorn temperature by using the microcanonical ensemble method, 
and investigated whether T = becomes stable near the Hagedorn temperature. 

When all the radii in the directions transverse to the Dp-Dp system are about the 
string scale, the result of the D + d = 9 case is in sharp contrast to that of the D + d < 8 
case. For D + d = 9 case, the sign of the coefficient of the |T| 2 term changes from negative 
to positive at slightly below the Hagedorn temperature. This implies that a phase transition 
occurs at this temperature and the Dp-Dp system becomes stable above this temperature. 
For D + d < 8 case, the coefficient remains negative near the Hagedorn temperature, so that 
such a phase transition does not occur. We thus concluded that the Dp-Dp pair, which is 
extended in all the non-compact directions, is created near the Hagedorn temperature. 

When some of radii in the directions transverse to the Dp-Dp system are much larger 
than the string scale, the result of the D + d = 9 case is also in sharp contrast to that of 
the D + d < 8 case. In D + d = 9 case, a phase transition occurs if E reaches (4.13) which 
satisfies (4.36). Therefore, a phase transition occurs for large enough E even if Rq is much 
larger than the string scale and the Dp-Dp system is not extended in all the (^-directions. 
Whether phase transition occurs or not depends on the topology of the target space. Strings 
are sensitive to the topology of the target space. In D + d < 8 case, there is no possibility 
that a phase transition occurs. 

Since our calculation is based on the work of Andreev and Oft [11], we must confirm 
the validity of their work. However, we only need to analyze the vicinity of T = in order to 
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investigate whether the phase transition occurs or not. Our calculation is valid if the mass 
square is shifted as equations (2.12) and (2.13) for small \T\. 

The conclusion of our previous work is that only the spacetime-filling D9-D9 pair 
is created near the Hagedorn temperature in the non-compact background [10]. Since the 
D9-D9 pair always satisfies D + d — 9, it is also created for large enough E in the toroidal 
background. Thus, speculation about the Hagedorn transition is not denied in the toroidal 
background case. 

The application of our model to the inflationary cosmology would be very interesting 
[12] [13]. If E is larger than (4.13) and satisfies the condition (4.36) in the early universe 
then the brane-antibrane pair is stable and the universe expands inflationary. But when 
E becomes smaller than (4.13) as the universe expands, the brane-antibrane pair becomes 
unstable and the tachyon starts to roll down from T = [24]. It must be noted that if we 
consider the dilaton gravity, the tension depends on time, because the tension of a D-brane 
depends on the dilaton field. 

Finally, the results that a brane-antibrane pair is created only when it is extended in 
all the non-compact directions is favorable to the Plank solid model of Schwarzschild black 
holes [25]. It might be interesting to study the relation between our model and the black 
holes. 
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